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Abstract. In this paper we present the solution to the problem of recovering rather arbitrary 
integral operator based on incomplete information with error. We apply the main result to obtain 
optimal methods of recovery and compute the optimal error for the solutions to certain integral 
equations as well as boundary and initial value problems for various PDE’s. 
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1. Introduction. Solutions to boundary (or initial) value problems for various 
partial differential equations require knowledge of a boundary (or initial) function. 
However, often time, those functions are not fully known and only partial information 
about them can be measured, e.g. values at some finite set of points, average values 
over small measurement intervals, values of N first consecutive Fourier coefficients, 
etc. Thus, it is very important to find an approximate solution based on available 
information on the boundary (or initial) function. Furthermore, it is also natural and 
important to develop methods that provide an optimal (in some sense) approximation 
to the true solution. These research questions have been explored under the theory 
of optimal recovery of functions and operators, which is an area of Approximation 
Theory that started to develop in 1970s. More information on the development of the 
area can be found, for instance, in [20, 28, 21, 12, 17, 24, 25, 11]. 

As for specific applications to recovering solutions of boundary and initial value 
problems, Magaril-Ill’yaev, Osipenko, and co-authors (see, for instance, [16, 22, 18]) 
have considered the problem of optimal L 2 -appi'Oximation of the solution to the Dirich- 
let problem for Laplace’s and Possion’s equations in simple domains (disk, ball, an¬ 
nulus) based on the first N consecutive Fourier coefficients of the boundary function 
(possibly given with an error). In order to solve this problem they have used methods 
of Harmonic Analysis and general results from Optimization Theory. 

In this paper we address related questions of optimal approximation of the solu¬ 
tion to several types of integral equations, boundary and initial value problems for 
PDE’s. We begin by solving a more general problem of recovering a rather arbitrary 
integral operator and sum of operators. We then present the optimal method of recov¬ 
ery as well as the optimal error. Next, we apply this general result to recover solutions 
to various boundary and initial value problems. Moreover, we present optimal meth¬ 
ods of recovery of the solution to boundary-value problems based on this incomplete 
information with error. Naturally, the solution to the problem when information with 
error is used will also lead to the solution to the problem with exact information. 
In this paper we focus on considering the Volterra’s and Fredholm’s linear integral 
equations as well as boundary value problems for wave, heat, and Poisson’s equa- 
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tions. Nevertheless, the developed method is more general and can be applied to 
other similar problems. 

The paper is organized as follows. Section 2 contains necessary definitions and 
notation as well as the formulation and solution of the main problem. In Section 
3, we solve the problem of optimal recovery of positive integral operators on classes 
of functions defined by moduli of continuity, based on information with an error 
about values of such functions at a fixed system of points. In Section 4, we use our 
general result from Section 3 to address optimal recovery problems for the solutions 
of Volterra and Fredholm integral equations of the second kind, systems of linear first 
order differential equations with constant coefficients, Poisson’s equation, the heat 
and wave equations. 

2. Statement and solution of the main problem. 

2.1. Definitions and notation. For l,m gN, we let be a collection 

of real linear spaces, be a collection of real linear normed spaces, and 

be a collection of real linear spaces. Set 

X := Xi X ... X Xm, y := Yi X ... X Yj, Z •.= Zi x ... x Zm- 

We write elements of spaces X, Y, Z as vector-columns, e.g. ir G X is a vector-column 
consisting of elements xi,X 2 ,..., Xm with xj G Xj, j = 1,... ,m. This allows us to 
equip spaces X, Y^ Z with natural coordinate-wise linear structure. In addition, in 
the space Y we introduce the norm 

(2.1) \\y\\y = \\y\U:=^{\\y,\\y,,...,\My,), 

where ip is an arbitrary norm in K*, monotone with respect to the natural partial 
order in R*. 

By 6 we denote zero of a linear space. It will be clear from the context what space 
is being discussed and, hence, we omit specifying it in the notation. 

Next, for a collection of linear operators Aij : Xj ^ Yi, i = 1,... ,1, and j = 
1,..., m, with domains of definition D (Aij) we consider operator matrix 



^Aii 

Ai2 ■ 



A 21 

A 22 

■ ^2m 

A := 

V^/i 

Ai2 

* -^Im J 


The matrix A defines the operator A : X ^ Y mapping an element x G X into 
the element y = Ax, which is a result of formal multiplication of matrix A by the 
vector-column x, i.e., for every i = 1,... ,1, the element yi is defined as 

m 

Vi — 

1=1 

For a given set of numbers aj G {—1,1}, J = 1, • ■ • j w, by a we denote the diagonal 
matrix 

fai 0 • • • 0 \ 

0 (72 • • • 0 

y 0 0 • ■ • am) 


a = diag((Ti,.. .,am) 
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We define the product of operator matrix A by matrix cr as a result of formal multi¬ 
plication of corresponding matrices, i.e. Aa is the operator matrix: 


Aa = 


/ CTi All 
(Tl A 21 


0'2Ai2 

172A22 


^mAirn \ 
A2m 


\(TiA/i (T 2 Aj 2 


^mAljji 




I 

In each of spaces Vj, j = 1,..., m, we select a class of elements Wj C n d{a,3), 

and consider the Cartesian products of classes Wj’s and their linear spans, respec¬ 
tively: 

W := Wi X ... X Wm C X, and spanW := spanWi x ... x spanIVm. 


Let us assume that a collection of operators Ij : spanW,- —>■ Zj, j = is 

given. We call information operators each of operators Ij and the operator matrix 
I = diag (/i,..., Im) ■ span W Z as well. Note that in majority of applications we 
consider in this paper, information operators will be linear. 

Finally, for a collection of operators Lj : Zj Xj, j = 1 , ..., m, by A L we denote 
the operator matrix obtained as a result of formal multiplication of operator matrices 
A and L = diag (Li, ..., Lm), where by AijLj we understand the composition of 
operators Aij and Lj. 

2.2. Optimal recovery problem and general lo-wer estimate for the error 
of recovery. In this paper we consider the problem of optimal recovery of operator 
A on the class W using information Ix on elements x from this class. 

For non-empty sets t/i,... Um in spaces Zi,..., Zm, respectively, we set 

U = Ui X ... X Um- 

We assume that instead of Ix, we know some element 1 G Ix + U, where U is given 
set, containing zero 0. In this situation we say that information Ix is known with 
U-error. Note that if U coincides with the origin of Z then ^ = Ix, and we say that 
information Ix is given exactly. 

An arbitrary mapping ^ . Z ^ Y is called a method of recovery. Given the 
operator A, class W, information / with [/-error, we define the error of recovery of 
operator A with the help of method $ as follows 

£■ (A; W;/;[/;$)= (A; VF;/;[/;$):= sup sup || Air — <i)z||^ , 

x^W 'z£lx-\-U 


and the error of optimal recovery of operator A as 

(2.2) S (A-,W;I;U) = S^(A;W;1;U) := inf _ £1 (A; W; I; [7; $) . 

The problem of optimal recovery of operator A: find the optimal error 
£ (A; W;I;U) and the method of recovery <!> (if any exists) delivering the inf in the 
right hand part of (2.2). 

Clearly, when U is the origin of Z, problem (2.2) reduces to the problem of optimal 
recovery of operator A on the class W based on exact information Ix on elements 

xgW. 
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Problems of optimal recovery of operators based on exact information were studied 
in [27, 5, 9, 20], and on approximate information in [15, 19, 21, 17, 16, 4]. We also refer 
the reader to the discussion of closely related questions in [28, 20, 12, 30, 24, 1, 2]. 

Let us provide the lower estimate on the error £ (^A;W; I;U). By W{U) we denote 
the following class: 

WiU) := {x£W : (Ix + U)n (I{-x) + U)^9} 

Proposition 1. If for every i = 1,..., / and j = 1,... ,m, the operator Aij is 
odd, W{U) ^ 0, and the class Wj is centrally symmetric, then 

(2.3) £4, {A;W-,1;U) ^ s_up_ px|| 

xeW(U) 


Proof. For any method of recovery $ : Z —>■ P, 
£,p {A;W;I;U;^) = sup \\Ax - ^z\\.ip 

r -sew 'i 
\ z€7x + U j 


> max ^ sup \\Ax — 4)2;||.0; sup || — Ax — ^z\\^ 
( T r -1 e w I 

. \ zeix+u J \ -e7(-x)+i7 / 


1 

^ 9 

^ xew(u)_^ 

^ sup 11^*1 
xeW(u) 


Ax — 4>z|| , + ll^a; + ^z\ 




ij) ' 


which completes the proof. □ 

In particular, when / is an even operator, the condition 

(7x + C7) n (T(-s) + U)^% 

is valid for every x W. Also, one can easily verify that there holds the following 

consequence from Proposition 1. 

Proposition 2. Let assumptions of Proposition 1 hold. In addition, for every 
j = 1,... ,m, we let Ij be odd, and Uj be centrally symmetric. Then 


(2.4) 


£.,p [A]W-, I;U) ^ sup ||Aa;| 

/ seW I 
\ Ix€U / 


■4>' 


Note that the lower estimate provided by Proposition 2 might not be sharp. 
Therefore, in the rest of this section we consider some general situations when in¬ 
equality (2.4) turns into equality. 

2.3. Normed lattices and positive operators. . Let us follow [26] in order 
to introduce the concepts of an ordered vector space, a normed lattice, and a positive 
operator. 

Definition 1. Given a linear space X over the field of real numbers R and a 
partial order “<x” on the set X, we call the pair {X, ^x) on ordered vector space if: 

1. X -<x y implies x + z -<x y + z, for all x,y,z € X; 

2. X -<x y implies Xx <x Xy, for all x,y € X and X G R+. 
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In what follows, for brevity (when it does not lead to confusion), we omit men¬ 
tioning partial order in the notation of an ordered vector space {X, ■<x)- 

Also, we reserve notation for the standard linear order in R. 

Definition 2. An ordered vector space X is called a vector lattice (the Riesz 
space), if any two elements x,y & X have supremum xV y := sup{a;; y} and infimum 
X Ay := inf{a;; y}. 

For a vector lattice X, by \x\ := a; V (—x) we define absolute value oi x & X. In 
addition, we call a norm on a vector lattice X a lattice norm, if \x\ -<x \y\ implies 
||x|| < ||y|| for all x,y € X. 

Definition 3. A normed lattice is a real normed space X endowed with an 
ordering “^x ” such that {X, -<x) Is a vector lattice and the norm on X is a lattice 
norm. 

Note that given a collection of ordered vector spaces Xj, j = their 

Cartesian product X is also an ordered vector space with respect to naturally defined 
partial order “-<y”: 

{x' (ij = 1 ,..., m, x'j -<Xj x'-). 

Similarly, for a collection of normed lattices Yi, i = 1,... ,1, their Cartesian product 

Y is also a normed lattice with respect to the norm || • ||^ defined by (2.1) and partial 
order 

Finally, we define the positive operator between ordered vector spaces as follows. 
Definition 4. Let X andY be ordered vector spaces. A linear operator T : X — >■ 

Y is called positive if 9 -Ay Tx whenever 6 -<x x. 

2.4. General results for positive operators. In this subsection we present 
some results on optimal recovery of positive operators and, in particular, identity 
operator. Furthermore, we show that under certain assumptions, once we know how 
to recover (in an optimal way) the identity operator on each of classes Wj, based on 
information Ij with 17,-error, we can recover (in an optimal way) any operator matrix 
A consisting of positive linear operators (and even operator matrix Aa) on the class 
W, based on information / with 17-error. 

Let X be a normed lattice, and by idx we denote the identity operator. We 
start with the problem of optimal recovery of the identity operator. Let Z be a real 
linear space, W C X he centrally symmetric class, / : X ^ Z be an odd information 
operator, and 17 C Z be non-empty centrally symmetric set. 

Proposition 3. If there exist an operator L : Z ^ X and a function p G W, 
lip G U, such that for any x £W and z G Z we have 

(2.5) {zgIx + U) => {—p^xx — Lz~<xip), 

then operator L is the optimal method of recovery of idx on the class W, based on 
information I with U-error, and 

(2.6) f {idx; W; I;U) = E {idx; W; I; U; L) = ||(^||^ . 

Proof. By assumption, for every x G W and z G Ix U we have 

—(f -<x X — Lz -<x p. 

Since X is a normed lattice, from the latter we obtain 


x-Lz\\x < \\p\\x- 
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Hence, 

£ {idx\W-J]U) ^ E {idx-,W\I\U]L) = sup ||x - ^ ||(^||x- 

x^W z£lx-\-U 

On the other hand lip € U. Hence, due to Proposition 2, we obtain 
£ {idx;W-,I-,U) ^ sup ||x||x > ||<p||x- 

r xew 1 

t ix€Crj 


□ 

Next, we present the result on optimal recovery of positive operators, which fol¬ 
lows from Proposition 3. 

Proposition 4. Under assumptions of Proposition 3, let Y he a normed lattice 
and A : X ^ Y he a positive linear operator with domain of definition n(H) D W. 
Then $ = AL is the optimal method of recovery of operator A on the class W, based 
on information I with U-error, and 

£ (H; W; I;U)=£ (H; W; /; U; $) = \\Ap\\y . 


Proof Indeed, let x € W and z & Ix -\- U he arbitrary. According to (2.5), we 
have —p -<x x — Lz -<x P- Hence, due to positivity of operator A, we obtain 

—Ap -<Y Ax — ALz -<Y Ap. 

Taking into account that || • ||y is a lattice norm, we deduce that 

\\Ax - ALz\\y ^ \\Ap\\y 

and, therefore, 


£ (A; W; P,U)^£ (A; W; /; U; AL) ^ || A^||y. 

The opposite inequality follows from Proposition 2. □ 

Remark 1. In Proposition 4 the condition that X is a normed lattice can he 
relaxed to the following one: X is an ordered vector space. 

Finally, we present the generalization of Proposition 4 to the case of optimal 
recovery of operator matrices. 

In order to state the corresponding result, we let Xi,... ,Xm be ordered vector 
spaces, Zi,..., Zm be real linear spaces, and Yi,... ,Yi be normed lattices, m,l gN. 
In addition, let A^ : Xj —>■ Yi, j = I,...,m and i = l,...,l, be positive linear 
operators with domains of definition V^Aij), and Ij : Xj —>■ Zj, j = l,...,m, be 

i 

odd information operators. Let also Wj C H 2? (Ay), j = l,...,m, be centrally 

i—1 

symmetric classes, and Uj C Zj, j = l,...,m, be centrally symmetric sets. Finally, 
let i/’ : K* —^ K be an arbitrary norm monotone with respect to the natural partial 
ordering in Rb 

Theorem 2.1. If for every j = 1,... ,m, there exist an operator Lj : Zj Xj 
and a function pj G Wj, IjPj G Uj, such that for any x G W and z G Z, we have 

{z g7x-\-U) ^ (Vj = 1,..., m, -pj -^Xj Xj - LjZj Pj) , 
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then for every a = diag{ai,... ,am), cFj G j = the operator 

$ = AaL is the optimal method of recovery of operator Aa on the class W based 
on information I with U-error, and, furthermore, 

{Aa;W-,T,U) = {Aa-W-l;U;<^>) = \\Alp\\ 


Proof. Let x &W and z & Ix -\- U. Then for any j = 1,... ,m, we have 

Xj - LjZj ipj. 

Since operators Aij are linear and positive, we obtain 

-^ijPj ^Yi -^ij(.Xj ^Yi 

From the latter we conclude that VtTj G {—1,1} 

^Yi ^j^ijXj ^Yi 

Summing up these inequalities over j, we see that 

m m m m 

— Aijipj -<Yi <XjAijXj — OjAijLjZj <Yi ^ijPj- 

1=1 1=1 1=1 1=1 

Therefore, 


CTl Aij Xj <Jj Aij L. 

1=1 1=1 


Zo 


~<Yi ^ijPj- 
1=1 


Since Il'Hy. is the lattice norm, from the latter we derive 


m m 


m 

ajAijXj — ajAijLjZj 


Hii (fj 

1=1 1=1 

Yi 

1=1 


and since $ = AaL, 


Hence, 


\Aax — HcrLzIl , = IlHcra; — <I>z|| , < ||Hw|| , . 

I W'l/j M 


£^{Aa-W-I-U-^) ^£j,{Aa-W-I-U) ^ \\Ap\\^. 


In order to obtain the lower estimate, we observe that aip G W and I{aip) G U. Due 
to Proposition 2 we obtain 

£jj {Aa;W;l;U) ^ sup ||Hcrs|| > ||Hcr(CT^)|| = ||H^|| 

/xew_\ ^ ^ 

\ IxGU J 


as a a = diag (1,..., 1) is the identity matrix. □ 
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3. Optimal recovery of integral operators. In this section we introduce the 
concept of an integral operator and apply Theorem 2.1 to the problem of optimal 
recovery of positive integral operators on classes of functions defined by moduli of 
continuity, based on information with an error about values of such functions at a 
fixed system of points. 

3.1. Integral operators on metric spaces. We follow [3] (see also [13]) to 
introduce the notion of integral operators on metric spaces. First, we let (M, be the 
space with cr-finite measure, i.e. M is some set and /r is a cr-finite measure on a-algebra 
Tjm of subsets in M. By 911 (M,/r) we denote the space of all /r-measurable jjL-a.e. 
finite functions defined on M (identifying //-equivalent functions as usual). The space 
911 (M, //) is equipped with the natural partial order : for every a;', x" G 911 {M, //), 

{x'< x") <t=> (for//— a. e. 1 G M, |a;'(l)| ^ |x"(l)|). 

Hence, we can consider the space 911 (M, /t) as an ordered vector space. 

Definition 5. (see [3], [13, Ch. 1, %2]) Let and {N,v) he spaces with 

a-finite positive measures, E be a linear manifold in A linear operator 

T : E ^ dJt{N, v) is called integral operator if there exists v x fi-measurable function 
K{s, f) such that for every x € E, 

(3.1) Tx{s) := / K{s,t) x{t) dp.{t), for pi —a. e. s G N. 

Jm 

The above integral is understood in the Lebesgue sense. The function K(s,t) is called 
the kernel of operator T. 

Remark 2. Clearly, an integral operator is positive if its kernel is v x p-a.e. 
non-negative. 

Next, we let M = Mp be a metric space endowed with the metric p. By Sp we 
denote the Borel ct- algebra of subsets of Mp, i.e. the minimal cr-algebra generated by 
open sets in Mp. We consider an arbitrary non-negative cr-finite measure p, : Ep ^ R+. 
For convenience, if the equivalence class in 911 {Mp, p) contains a continuous function, 
then we identify this function with the whole equivalence class. 

By Bp and Cp let us also denote the sets of //-essentially bounded and p-a.e. 
continuous functions x : Mp —>■ K, respectively. 

For a compact set M' C Mp, we let xm' stand for the characteristic (or indicator) 
function of the set M'. We consider classes 

Bp{M') := {x G Bp : suppx C M'} , 

Cp{M') := {x = y xm' ■ V G Cp] , 

Cp{M') := {x G Cp : suppx C M'} . 

By definition, Cp{M') C Cp{M') C Bp{M'). 

3.2. Classes and generalized Voronoi cells. We recall that a function 
uj : K+ —K+, R+ := [0, oo), is called a modulus of continuity (see, for example, [10]) 
if iu;(0) = 0, a; is continuous, non-decreasing, and semi-additive function. The latter 
means that uj{t' 1") ^ w(t') -I- w(t"), for every t', t" G M+. 

We consider the problem of optimal recovery of positive integral operators on the 
classes defined by a modulus of continuity w: 


H;(M') := {a; G Cp{M') : \x{t') - x{t'')\ ^ tv {p{t',t")), ^t',t" G M'} , 
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In addition, we assume that information on functions x € H‘^ {^') or x € H‘^ {M') 
is provided by information operators I : {M') —>• M", n S N, of the form 

Ix = Iqx := {x {qi ),.. .,x{qn )), x € C^[M'), 

where Q = {qj}^^i is a hxed set of points in M', and is known with Ug-error, e € R", 
where 


Ug — [ Cl, Cl] X ... X [ Cyi]. 


Next, we construct the operator L that would satisfy assumptions of Proposi¬ 
tion 3. To this end, we first introduce the following two functions: 


(3.2) 


and 


T{t) = T^,Q,g{t) 


min {Cj + Oj{p{t, qj ))), t G M', 

j — l,...,rL 

0, t€M\M', 


(3.3) f(t) = fu,,Q,git) := min{Tu,,Q,g{t); uj{p{t,dM'))} , t e M, 


where p{t, dM') := inf p(t, s) denotes the distance between point t G M and the 

seSM' 

boundary dM' of M'. One can easily verify that r G f G and 

Iqt = Iqt G Ug. 

Next, we define generalized Voronoi cells. To this end, we first let 
fio = no(w, Q, e) := {t G M' : uj{p{t, dM')) ^ T^,Q,g{t)} , 


and, for j = 1,..., n, we consider 

n' = n'(w, Q, e) := {t G M' : r^.Q,e(<) = ej + w{p (t, g^))} . 

Then, generalized Voronoi cells (see Figure 3.2) are defined iteratively as follows 

j-i 

ni:=ni, nj:=n'\(Jn^, j = 2,...,n, 

and 

ni = nj\no, j = i,...,n. 

Since every function x G C^(M') is ^-measurable on M in the Borel tr-algebra Ep, 
we see that sets Hi, 112,..., n„, and sets Ho, Hi,..., 11,1 are /i-measurable. Moreover, 
due to construction, we conclude that both collections of sets are pairwise disjoint 
and 


M' = Hi u... u n„ = fio u ill u... u n„. 

Finally, we define operators L = : R" —>• Bp (-^0 and L = : R" 

Bp (M') (see Figure 3.2) as follows: for every z G M", 

/ Zj, t G Hj, j = 1,... ,n, 

\ 0 , teM\M', 


Lz{t) : 
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Fig. 1. Generalized Voronoi cells on the segment [a^b] 


+ *3 



Fig. 2. Method of recovery of a function defined on the segment based on information about its 
values at 6 points with error 


and 


t G n^, j = 
t€UoU{M\M'). 

Note that when some of Cj ’s are large enough, it is possible that some of the sets 
Ilj and Ilj are empty. This, in turn, means that information at the corresponding 
point Qj is “ignored” by operators L and L. 

3.3. Optimal recovery of operators on the class . In this section, we 
present some important direct corollaries from Theorem 2.1 for positive integral op¬ 
erators and their sums. 

For a space (M, fi) with cr-hnite measure, we let Li (M, be the space of abso¬ 
lutely integrable functions x : M M. with the standard norm 

\\x\\i= [ \x{t)\dn{t). 

J M 

We start with the corollary from Proposition 4 for positive integral operators. In 
order to state the corresponding result, we let n G N, w be a modulus of continuity, 
e G R", Mp be a metric space, ^ be a ct- finite measure on Sp, M' be a compact 
subset of Mp, (iV, v) be a space with cr-finite measure, A : Bp{M') —> 911 {N, v) be an 
integral operator with the non-negative kernel K, Y C DJI {N, v) be a normed lattice 
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such that A {M')) CY, Q = C M' be a fixed system of points. 

Theorem 3.1. Let either W = ip = T^^Q^e> L = L^ Q f^, or W = 

p = L = Then the method = AL is the optimal method of 

recovery of operator A on the class W based on information Iq with Ue- error, and 


(3.4) f (A; W; Iq; U^) = £ (A; W; lQ;Ue; <&) 



K(-, x) p(x) dp(x) 


Y 


Proof We consider only the first case when W = p = Ta,,Q,e, L = Luj,Q,e 

as the proof of the second case follows similar arguments. 

We let X G and z € Ix + U he. given, and let IIi,..., n„ be generalized 

Voronoi cells on the set M'. We observe that for every j = 1,..., n and t G 11^, 

x(t) - Lz{t) = \x{t) -Zj\^ \x{t) - x{qj)\ + \x{qj) - Zj\ 

< qj)) + ej = Tuj,Q,e{t) = p(t). 

Hence, we conclude that —p -< x — Lz -< p where is the natural partial order 
in dJt {Mp, p). Since H is a linear positive operator, conditions of Theorem 2.1 are 
satisfied, and by applying it, we complete the proof. □ 

When Y = Li {N, iz), we obtain the following consequence from Theorem 3.1. 
Corollary 1. Under assumptions of Theorem 3.1, we take kernel K to he 
V X p-integrahle and Y = Li{N, v). Then 

(3.5) £ {A;W;lQ;Ue) = [ p{x) [ K{y,x)dv{y)dp{x). 

Jm' J n 

In particular, when the inner integral in (3.5) is independent of x and is denoted by 
Ck, we have 


S {A;W]lQ;Ue) = Ck [ 

JM' 


Next, we state the consequence of Theorem 2.1 for the problem of optimal recovery 
of sums of positive operators on classes defined by moduli of continuity. We need the 
following notation. Let m,l G N and i/) : K* —>■ R be a norm monotone with respect to 
the natural partial order in M.K For every j = 1,..., m, we let Uj G N; ujj be a modulus 
of continuity; Cj G Rp’ and Uj = ; Mj = Mp^ be a metric space, pj be a cr-finite 

measure on Yp .; Mj be a compact subset in Mj ; Qj C Mj be a fixed systems of nj 
points, and Ij := . For every i = we let (W, t'i) be a space with cr-finite 

measure. For every i = 1,... ,l and j = 1,..., m, we let Aij : Bp. (M'j) —> 911 (W, Ui) 
be an integral operator with non-negative kernel Kij, and Yi C DJt (W, Vi) be a normed 

m 

lattice such that IJ Aij [Bp. {M'j)) C Yi. 
i=i 

Theorem 3.2. Let, for every j = either Wj = Hp) {M'j), Pj = 

Tajj,Qj,ej, Lj = or Wj = Hp) {Mj), Pj = Lj = Lt^.^Q.^g.. Then, 

for every a = diag (cti, ..., am), where aj G {—1,1}, j = 1, ..., m, the method $ = 
AaL, where L = diag {Li ,..., Lm), is the optimal method of recovery of operator Aa 
on the class W based on information I with U-error, and 

£^(Aa;W;l;U)=£^(Aa;W;l;U;^) = . 


(3.6) 
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The next proposition deals with the optimal recovery problem of the sum of 
integral operators in the space Li{N, v). 

Corollary 2. Let assumptions of Theorem 3.2 hold. In addition, we assume 
that kernels Kij are Vi x fij-integrable (i = 1,... ,1, and j = \,... ,m),Yi = Li (Ni, Vi), 
and if is £i-norm on Then 

(Aa; W; T; U) = (Aa; W; T; U; $) 

I m p p 

= V V / / {y, x) (fij (x) dyj{x) dvi{y). 

i=i j=i Jm', 


4. Applications. In this section we demonstrate how main results of this pa¬ 
per can be applied to the problems of optimal recovery of the solutions to integral 
equations, and boundary and initial value problems for differential equations. 

In Section 4.1 and 4.2 we present optimal methods and error of recovery of solu¬ 
tions for Volterra and Fredholm integral equations of the second kind. 

We then present optimal methods and errors of recovery for solutions of systems of 
linear first order differential equations with constant coefficients, Poisson’s equation, 
the heat and wave equations. Certainly, our approach is not restricted to optimal 
recovery of the solutions to mentioned equations and is applicable to a wider range of 
integral equations, ODE’s, and PDE’s. 

4.1. Linear Volterra integral equations of the second kind. Let [a, &] C R 

be the interval, y be the Lebesgue measure on [a,b], function / : [a,b] ^ R and 
kernel k : [a, b]^ —R be given, x : [a, ^ R be unknown function. A linear Volterra 

equation of the second kind is the equation 

(4.1) x(t) = f{t) + f k{t, s) x{s) dy{s), t€[a,b]. 

J a 

We let P be the resolvent kernel for k: 

CO 

r(t, s) := y] s), t,s€[a,b], 

n—1 


where fci = k, and, for n = 2, 3,..., 

kn{t,s) := / k{t,T)kn-i{T,s) dy{T), t,s€[a,b]. 

J S 

It is well known (see, for instance, [14, Theorem 3.3]) that for continuous function 
/ and kernel k the solution to (4.1) exists, is unique, and can be written in the form 

(4.2) x{t) = f{t)+f T{t,s) f{s)dy{s), t€[a,b]. 

J a 

Let w be the modulus of continuity, Q be the set of n € N points on [a,b], 
and e G R". Let us consider the problem of optimal recovery of the solution to 
equation (4.1) under assumptions that the values of function / G Hff ([a, &]) are known 
at the system of points Q with t/g-error. 
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By (4.2), the solution x to problem (4.1) can be considered as the image x = Af 
of the function / under the operator which is the sum of identity operator and an 
ntegral operator with the kernel K{t, s) = k{t, s) ■ X{a,t){s)- 

Let Y C 91t([a, &],/r) be a normed lattice containing the set of i3^([a, 6]) of fj,- 
essentially bounded on [a, b] functions. Then in view of Proposition 4, there holds 
true the following 

Corollary 3. If the kernel k is non-negative on [a, 6]^ then the method $ = 
ALi^^Q^e is the optimal method of recovery of operator A on the class H‘^([a,b]) based 
on information Iq with Ue- error, and, furthermore, 

£{A-H^{[a,b]y,lQ-,Ue) = £{A-H^{[a,b]y,lQ-,Ue;<I>) = \\AT^,Q,e\\Y 

rk) 

= t^,qA-)+ r((-),s)T<^,Q,e(s)d/r(s) 


It follows directly from Corollary 3 that for Y = Li{[a, b],p,), 


£{A-HA[a,b]yiQ-,Ue) 


T{t, s) dpL{t) 


Tu;,Q,e{s) 


4.2. Linear Fredholm integral equations of the second kind. Let [a, b] C 

M, /r be the Lebesgue measure on [a,b], / : [a, 6] —K be continuous function, x : 
[a, ^ M be unknown function, and kernel fc : [a, 6]^ —>■ R be such that 


(4.3) 


\k{t,s)\^ dp,{t) dAs) < 1- 


The linear Fredholm integral equation of the second kind is the equation 


(4.4) 


x{t) = f{t)-\- f k{t, s) x{s) dAs)- 


By r we denote the resolvent kernel for fc: 


OO 

r(t, s) := ^ fc„(t, s), t,s€[a,b], 

n—1 


where fci = k, and, for n = 2, 3,..., 

kn{t,s)= / kn-i{t,u) k{u, s) dAAi t,sG[a,b]. 

J a 

It is well known (see, for instance, [6, p. 44]) that the unique solution to (4.4) is given 
by 

(4.5) x[t) = f{t)+[ T{t,s) f[s)dAs), t€[a,h]. 

J a 

We let w be a modulus of continuity, Q be a set of n S N points on [a, b], e G R" . 
We consider the problem of optimal recovery of the solution to equation (4.4) under 
assumption that the values of function / G &]) are known at the system of 

points Q with [/g-error. 
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In virtue of (4.5), the solution x to problem (4.10) can be considered as the image 
X = Af of the function / under the operator A, which is the sum of identity operator 
and an integral operator with the kernel K{t, s) = kit, s). 

Let Y C 91l([a, &], /r) be a normed lattice containing the space B^{[a, 6]). Then by 
Proposition 4, there holds true the following 

Corollary 4. If the kernel k is ^-a.e. non-negative on [a, b]^ and satisfies (4-3), 
then the method ^ = A Lu>,Q,e is the optimal method of recovery of operator A on the 
elass H‘^{[a,b]) based on information Iq with Ue-error, and, furthermore. 


f (A; iJ-([«,&]);J q; C/e) = £ {A-H‘^{[a,b])-lQ-,Ue;‘^) = \\At^,q,, 


Tuj,qA-)+ / r((-),s)r,^,Q.e(s)dct(s) 

J a 


Y 


It follows directly from Corollary 4 that for Y = Li([a,6],/i) 


£ {A-HA[aA)AQ-Ye) 


T{t, s) dpL{t) 


Tui,Q,e{s) dfl{s)- 


4.3. Optimal recovery of solutions to the systems of differential equa¬ 
tions. Next, we consider the problem of optimal recovery of the solution to initial 
value problem for the system of linear first order differential equations with constant 
coefficients. Let us introduce several notation: let d S N, S' be d x d matrix with 
real entries, [a, 6] C M be a finite interval, p, be the Lebesgue measure on [a, b\. By x 
we denote vector function x = {xi ,..., Xd) consisting of functions Xi GiXft ([a, &], /i), 
i = 1,..., d. Finally, we let g be a continuous function, and p € be some point. 

The system of linear nonhomogeneous equations has the form: 


f x'(t) = S x(t)-h q(t), tG[a,b], 
\ x{a) = p. 


It is well known that the solution to (4.6) exists, is unique, and is provided by: 

(4.7) x{t) = P + [ q{u) d^u), t S [a, b], 

J a 

where e^ stands for the exponent of matrix M which is the series 


e" := y: 

3=0 


j! 


Next, we let Wi,..., be given moduli of continuity, Wi := and 


W2 X ... X id-''([a, 6]), 

Q be the given set of n G N points on [a,b], Ii := id^d : —)■ R'^ and I 2 ■= 

diag {Iq, ... ,Iq), consisting of d operators Iq : C'^([a, 6]) ^ R", be information 
operators, e G R^J. and ei,..., G R" be the errors describing information, Ui := Ue 
and U 2 := x ... x Ug^^- 

Let us consider the problem of optimal recovery of the solution to the system (4.6) 
under assumptions that the initial value p is known with C/i-error, and, for i = I, ..., d. 
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the values of component Qi of the function q £ W 2 at the system of points Q are known 
with C/ei-error. 

Next, for by t £ [a,6] we denote the element of matrix 

^S{t-a) located in the tth row and the jth column. In addition, we consider operators 
Bij : R —>• C^([a, &]) and &]) —>■ &]) defined as follows 

(Bijx) (t) = kijit) ■ X, t £ [a, &], X £ R, 

and 

{Cijx){t)= / kij{u) x{u) dfi{u), t€[a,b], x G Bfj_{[a,b]). 

J a 

In view of (4.7), the solution x to (4.6) is the sum x = Bp+Cq of images of initial value 
p and function q under operator matrices B and C, respectively. This observation 
allows us to consider the problem of optimal recovery of the solution to equation (4.6) 
as the problem of optimal recovery of the matrix operator A = {B C) on the class 
W = Wi X 11^2 based on information / = diag(/i,/ 2 ) with U = Ui x C/ 2 -error. 

Next, we remark that a square matrix S is called essentially non-negative if every 
non-diagonal entry of this matrix is non-negative. It is well known that for such 
matrix S', operators Bij and Cij (the entries of matrix operator A) are positive. 

For i = I,..., d, we let Yi C dJl{[a,b], p) be a normed lattice, containing space 
i3^([a, 6]), and let ip he a monotone norm in R'^. Let also Li := idg^d, and L 2 ■= 
diag (La;i,Q,ei, • ■ ■, LuJd,Q,ed)- III addition, we set ri := e and 

12 ■ ('fivi,Q,ei I ■ ■ ■ ! 

Applying Theorem 2.1, we obtain the following 

Corollary 5. Let S be essentially non-negative dx d matrix. Then ^ = AL, 
L = diag{Li,L 2 ), is the optimal method of recovery of operator A on the class W 
based on information I with U-error, and, furthermore, 

{A;W;T,U) = {A;W-,T,U;^) = p t||^ = ||;Be + Ct2||^ , t = ti xts. 

One can easily adjust the arguments of this section to solve the problem of optimal 
recovery of the solutions to the system of linear homogeneous equations and to the 
system of linear nonhomogeneous equations with homogeneous initial values. 

4.4. Optimal recovery of the solution to the Dirichlet problem for Pois¬ 
son’s equation. In this section, we consider the problem of optimal recovery of the 
solution to the Dirichlet problem for Poisson’s equation. We let d £ N, | • | denote the 
standard norm in Euclidean space p he the standard Lebesgue measure in R'^, 
D C R'^ be a bounded domain with C^-boundary, a be the surface area measure on 
the boundary d£l. As usual, A stands for the Laplace operator. In addition, we let 
C{£t U dil) be the space of continuous on D U dil real-valued functions, and 
be the space of continuous functions having continuous first, and second order partial 
derivatives inside LI. 

The Dirichlet problem for Poisson’s equation consists of finding a function x £ 
C{Ll U dLl) n (7^(17), called the solution, which satisfies 


(4.8) 


-Ax{t) = f{t), t £ D, 
x{t) = g{t), t £ dLl. 
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It is well known (see [7]) that the solution x to (4.8) exists, is unique, and can be 
presented in the form 

(4.9) x{t) = [ G{t,s) f{s)d^i{s) - f ^{t,s)g{s)da{s), t S O, 

Jn Jan on 

where G{t,s) is Green’s function of the domain fl, and is the outer normal 

derivative of G. 

Below, we assume that the domain and it boundary dfl are endowed with 
the respective Euclidean metric, and the metric which agrees with the surface area 
measure a on dil. We consider classes Wi := and W 2 '■= i7“^(9f2), where 

a;i,a ;2 are given moduli of continuity, finite sets of points Qi C and Q 2 C d^l 
consisting of, respectively, ni € N and n 2 € N points. We also assume that information 
operators Ii := Iq^ and I 2 '■= IQ 2 , and sets Ui := Ue^ and U 2 ■= Ue 2 , where ei G 
and 62 G R"^, describing the error of information, are given. 

Let us consider the problem of optimal recovery of the solution to the prob¬ 
lem (4.8) under assumptions that the values of functions / G Wi and g G W 2 at 
systems of points Qi and Q 2 are known, respectively, with Ui and [/ 2 -errors. 

By (4.9), the solution x to (4.8) is the sum x = Aif + A 25 of images of functions 
/ and g under integral operators Ai : -G VJl {ilU dft, g) and A 2 : Bcr(dfi) -G 

(fl U dil, /i) respectively with kernels 

Ki(t, s) = G(t, s), t,sGfl, 

and 

dG 

K2{t,s) = ten, s G dfi. 

Hence, the problem of optimal recovery of the solution to the problem (4.8) can be 
reformulated as the problem of optimal recovery of the matrix operator A = {Ai A 2 ) 
on the class W = WiX W 2 based on information / = diag(/i,/ 2 ) with U = Ui x U 2 - 
error. 

Since both operators Ai and A 2 are positive, the assumptions of Theorem 3.2 are 
satisfied. For convenience, we let F C SOI (H U Oft, g) be a normed lattice, containing 
space B ^ (flUttn), , and z — 1;2. 

Corollary 6. The operator ^ = AL, with L = diag {Li, L 2 ), is the optimal 
method of recovery of operator A on the elass W based on information I with U- 
error. Moreover, the optimal error is 

£{A;W-J-U) = £ {A;W-J;U-^) = Pt||.^, r = n x r2. 


It follows directly from Corollary 6, that for Y = Li (H U d£t, g) 

£ [A;W; I;U) = j ri(s) j G{t,s) dg{t) dg{s) 

Jq j q 

- [ '^2(5) [ -j—{t,s)dg{s)da{s). 

Jan Jn on 

In particular, when is the disk of radius r centered at the point a G R"^, we have 

£ {A;W-J;U) = ^ [ {r'^ - \s - a\'^) ti{s) dg{s) + ^ [ T2{s)da{s). 

^ Jn ^ Jan 


Optimal recovery of integral operators and its applications 


17 


One can apply similar arguments to solve the problem of optimal recovery of 
solutions to the Dirichlet problem for Laplace’s equation, and to the homogeneous 
Dirichlet problem for Poisson’s equation. 



Fig. 3. Solution to the Laplace equation 



Fig. 4. Recovery of the solution to the Laplace equation 


4.5. Optimal recovery of the solution to the initial value problem for 
heat equation. The heat (or diffusion) equation describes the evolution in time of 
the density of some quantity such as heat, chemical concentration, etc. In the present 
section we apply results on optimal recovery of integral operators to the problem of 
optimal recovery of the solution to initial value problems for this type of PDE. 

In what follows, we use the following notation. Let D = x (0,+oo), d G N. 
Let also | • | stand for the Euclidean norm in M'’*. For a function a: : U —>■ M, we denote 
by Ax its space-coordinates Laplace operator, i.e. 



In addition, we let V to be a manifold in D. By we denote the standard Lebesgue 
measure corresponding to the dimension of the manifold, equipped with this measure, 
i.e. for the space SOI (H,/r) measure ^ is (d + l)-dimensional Lebesgue measure, 
while for the space 9Jl the measure fi is d-dimensional Lebesgue measure. We 

take / : H > R to be a continuous, compactly supported function with continuous 
partial derivatives and , j,k = 1,..., d. Finally, let g : R.'^ ^ R be a 

continuous, bounded function. 









18 


Babenko V., Babenko Y., Parfinovych, Skorokhodov 


The Cauchy problem for the heat equation consists of finding a function x € 
C{D U dD) n C‘^{D), called the solution, which satisfies 


\ x{v,0) = g{v), V 


It is well known (see, for instance, [ 8 ]) that the solution x to (4.10), satisfying a growth 
condition 

(4.11) |a;(u, t)| < , {u,t) € D, 


for some constants a,/? > 0 , exists, is unique, and can be presented in the form 

e^p(-feS) 


(4.12) 


x{u,t) = 


f! 

Jo Jr“ 


/o JRrf ( 47 r(t - s ))'^/2 

»p(-‘= 9 *) 


/R'i 


(47rt)'^/^ 


f{v,s)dg{v,s) 
g{v)dg,{v), {u,t) e D. 


Below, we assume that D and are endowed with the Euclidean distance. Mi C 
D and M 2 C are compact sets. We consider the classes Wi := and 

W 2 ■= M“^(M 2 ), where uji,oj 2 are given moduli of continuity, finite sets of points 
Qi C Ml and Q 2 C M 2 consisting of, respectively, ni G N and n 2 G N points. We 
also assume that information operators Ii := Iq^ and I 2 '■= IQ 2 , and sets Ui := Uei 
and U 2 '■= Ue 2 describing the error of information, where ei G and 62 G R"^, are 
given. 

Let us consider the problem of optimal recovery of the solution to the prob¬ 
lem (4.10) under assumptions that the values of functions / G Wi and g G W 2 at 
systems of points Qi and Q 2 are known respectively with Ui and [/ 2 -errors. 

Due to (4.12), we can easily see that the solution x to (4.10), satisfying the 
growth condition (4.11), is the sum x = Aif + A 2 g of images of functions / and g 
under integral operators Ai : B^{D) 91t(D,/r) and A 2 : i3^(R'^) with 

respective kernels 


Ki{{u,t), {v,s)) 


exD ( 

exp J ^ 

(Mt - s))d/^ •X(o.t)(s 


{u,t), {v,s) G D, 


and 


K2{{u,t),v) 


=p(-^) 

(4M)'^/2 ’ 


{u, t) £ D, V £ R"^. 


Hence, the problem of optimal recovery of the solution to problem (4.10) satisfying 
growth condition (4.11) can be reformulated as the problem of optimal recovery of 
the matrix operator A = {Ai A 2 ) on the class W = Wi x W 2 based on information 
I = diag ih.h) with U = Ui x [/ 2 -error. 

Since both operators Ai and A 2 are positive, the assumptions of Theorem 3.2 
are satisfied. For convenience, we let Y C 31l(iV, g) be a normed lattice containing 
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the space of /i-essentially bounded and integrable on N functions, Li = and 

'^U)i,Qi,ei ; ^ f ; 2. 

Corollary 7. The operator ^ = AL, L = diag{Li,L 2 ), is the optimal method 
of recovery of operator A on the class W based on information I with U-error. More¬ 
over, 


£ {A-W-I-U) = £ {A-W-I-U-<^>) = llArlly, r = Ti xt2. 

In particular, for Y = Li{N,fx) we obtain 

1. If iV = X {to} where to > 0, then 

£ {A]W;7;U) = f [ ti{v,s) dp{s)dp{v) -\- f T 2 (n) d/i(n); 
Jr** Jo Jr-^ 

2. If TV = {uq} X (0, +oo) where uq € is fixed, then 


£{A-W-I-,U) = 


4 d-ir(d _ 


-d/2 



T 2 (u) 


/R'i |wo — V 


R'i |mo — V 

d-2 Mv) >; 


- j^dp{v) dp[s) 


3. If TV = {moI X {to} where uq € R'^ and to > 0 are fixed, then 

£(A-,W-,T,U) = [ 

Jo Js. 


Ao r exp(-^2^)ri(?;,s) 


lo (47r (to - s))^/2 


dp{v) dp{s) 


exp - 




/r‘^ (iTrto)'^'^^ 


T 2 (i;) dp{v). 


Finally, we note that one can adjust the above arguments to solve the problem 
of optimal recovery of the solutions to the Cauchy problem for the homogeneous heat 
equation, and to the homogeneous Cauchy problem for the heat equation satisfying 
growth condition (4.11). 

4.6. Optimal recovery of the solution to initial value problem for the 
wave equation. In this section we consider the problem of optimal recovery of the 
solution to the wave equation. We recall that the wave equation describes wave 
propagation in a media and is a simplified model for a vibrating string (d = 1), 
membrane (d = 2), or elastic solid (d = 3). 

For the purpose of this section, we take d = 1,2,3. Let | • | be the Euclidean norm 
in R"^. We introduce notation of a d-dimensional ball Bt{u) C R"^ and a d-dimensional 
sphere St{u) C R'^”^ centered at point m G R"^ with radius t > 0. In addition, we let 
_D = R"^ X (0, -boo). Similarly to the previous section, for a function x : U —)• R, we 
denote by Ax its space-coordinates Laplace operator, let A C U be a manifold, and 
pL be the standard Lebesgue measure corresponding to the dimension of the manifold, 
equipped with the measure. Let also / : D —>■ R, g : R'^ —)• R, and : R'^ —> R be 
some functions. 

First, we let d = 1 and consider the Cauchy problem for one-dimensional wave 
equation, which consists of finding a twice continuously differentiable function x : 
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£)—>•]&, called a solution, satisfying the system of equations: 


(4.13) 


Xtt{u,t) - Xuu{u,t) 

x{u,0) =g{u), 
xtiu,0) = h{u), 


f{u,t), {u,t)eD, 
M e K, 

M S M. 


If / is continuous, g is twice continuously differentiable, and h is continuously differ¬ 
entiable, then the unique solution to the problem (4.13) is delivered by the Dalambert 
formula: 


(4.14) 


= - 


t pU-\-S 

f{v, t — s) dfi{v) dg.{s) 
g{u-t) + g{u + t) 1 


U-\-t 


h{v) dfi{v). 


As we are interested in applications of Theorem 3.2, we would need to assume that g 
and h are compactly supported and have a majorant for their modulus of continuity. 
This means that g and h are continuous, but might be non-differentiable. Hence, 
the formula (4.14) does not deliver the solution to the problem (4.13). Therefore, we 
follow [23] to introduce the concept of the generalized solution to the problem (4.13). 

Let g,h G (7^(0), where H is a compact in R. Let also be a sequence of 

twice continuously differentiable functions converging uniformly on H to the function 
g, as n ^ oo. Similarly, let be a sequence of continuously differentiable 

functions converging uniformly on to /i, as n —^ oo. The function x(u,t) is called 

the generalized solution to the problem (4.13) if it is a limit of uniformly converging 

sequence of solutions {a:("'^(u,t)}^^ to the wave equation x[^\u,t) — x^{u,t) = 
f(u,t), {u,t) S H, with initial conditions 

x^"^u,0) = g^’^^(u), u € R, 

Xf’^\u,0) = u e R. 


Using the above definition, the Dalambert formula (4.14) delivers the unique gener¬ 
alized solution to the problem (4.13). 

Next we let Mi C D, and M 2 , M 3 C R be compact sets. For i = 1,2,3, we let 
uJi be modulus of continuity, Wi := be a class of functions, Qi be the set 

of Ui G N points on Mi, li := /q. be information operator, and Usi be the error of 
information, e; € R"*. 

We see that the generalized solution x to (4.13) can be presented as the sum 
X = A\f + A 2 g + A^h of images of functions /, g, and h under linear positive operators 
Ai : H^(Mi) ^ A 2 : B^{M 2 ) ^ and A 3 : ^ 

This allows us to consider the problem of optimal recovery of the generalized solution 
X to the wave equation (4.13) as the problem of optimal recovery of operator A = 
(Ai A 2 A 3 ) on the class W = Wi x W 2 x W 3 based on information / = diag {Ii, 12 , 13 ) 
with U = Ui X U 2 X [/ 3 -error. 

Let us formulate the following corollary from Theorem 3.1. Let Y C dJl{N,g) 
be the normed lattice containing essentially bounded and integrable on N functions, 
and, for i = 1,2,3, we let Li := 

Corollary 8. The operator ^ = AL, L = diag {Li, L 2 , L 3 ), is the optimal 
method of recovery of operator A on the class W based on information I with U- 
error. Moreover, 

£ {A;W;7;U) = f (A; W; T; U; 4>) = ||Ar 


r = n X r 2 X T 3 . 
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In particular, for Y = Li{N,iy), and iV = M x {to}j where to > 0 is fixed, we 
obtain 


r^o p+oo 

£{A;W;T,U) = s Ti{v,to — s) d^j.{v) d^{s) + I T 2 {v)d^{v) 

Jo J — oo 

/ oo 

T 3 {v)dfJ.{v). 

-OO 


Next, we let d = 2,3, and consider the Cauchy problem for the wave equation with 
zero initial form, i.e. g = 0. This problem consists of finding a twice continuously 
differentiable function x : D ^ W, called a solution, which satisfies the system of 
equations: 

{ xtt{u, t) - Ax{u, t) = f{u, t), {u, t) G D, 

x{u, 0) = 0, u G R'^, 

Xt{u, 0) = h(u), u G R"^. 


If / is continuous, and h is twice continuously differentiable then the unique solution 
to the problem (4.13) is delivered by the Poisson formula when d = 2: 


(4.16) 


x{u,t) = - 


s) 

B,{u) 

f h(v) dg,(v) 


: dg,(v) dg,{s) 


27 ^ JBtiu) -\u-' 


and by the Kirchhoff formula for d = 3: 


(4.17) 


x{u,t) = ^ 


Bt(u) 

I 

4TTt 


f{v,t- \u-v\ 
|m — uj 

/ h{v)da{v), 
'St(u) 


■ dg{v) 


where a stands for the surface area measure of the sphere St (u). 

Similarly to the case d = 1, we would need to assume that h is compactly sup¬ 
ported and has a majorant for the modulus of continuity. Therefore, we follow [23], 
and introduce the generalized solution to the problem (4.15) as follows. 

Let h G where is a compact in R, and }„-i be the sequence of twice 

continuously differentiable functions converging uniformly on 11 to /i, as n —>■ oo. The 
function x(u,t) is called the generalized solution to the problem (4.15) if it is a limit 
of uniformly converging sequence of solutions (it, to the wave equation 

x[t\u,t) - X^uu {u,t) = f(u,t), {u,t) G D, with initial conditions 


x("^(it,0) = 0, It G R, 
x[^\u,0) = It G R. 


Using the above definition, the Poisson and the Kirchhoff formulas deliver the unique 
generalized solution to problem (4.15). 

Now, we let Mi C D and M 2 C R'^ be compact sets. For i = 1, 2, we let uji to 
be a modulus of continuity, Wi := be a class of functions, Qi be the set 

of rii G N points on Mi, li := Iq. be information operator, and Ue^ be the error of 
information, Ci G R"g 
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We see that the generalized solution x to (4.15) can be presented as the sum 
X = Aif + A 2 h of images of functions / and h under linear positive operators Ai : 

—>• dyi{N,^), and A 2 : —)• This allows us to consider the 

problem of optimal recovery of the generalized solution x to the wave equation (4.15) 
as the problem of optimal recovery of operator matrix A = {Ai A 2 ) on the class 
W = Wi X W 2 based on information I = diag (Iijh) with U = Ui x [/ 2 -error. 

From Theorem 3.2 we obtain the following corollary. Let Y C dyi{N, fi) be the 
normed lattice containing /r-essentially bounded and integrable on N functions, and, 
for z = 1, 2 , we let L* := n 

Corollary 9. The operator ^ = AL, L = diag{Li,L 2 ), is the optimal method 
of recovery of operator A on the class W based on information I with U-error. More¬ 
over, 


8 {A-W-I-U) = £ {A-W-I-U-^) = ||/It||^, r = Ti xr 2 . 

In particular, for Y = Li{N,p), and N = x {to}, where to > 0 is fixed, we 
obtain 

£ {A;W;I;U) = [ s f Ti{v,to - s) dp{v) dp.{s)-\-to f T 2 (v)dpL{v). 

Jo Jr'^ 
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